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The dynamics of collapsing and expanding cylindrically symmetric gravitational and matter fields 
with lightlike wave-fronts is studied in General Relativity, using the Barrabes-Israel method. As 
an application of the general formulae developed, the collapse of a matter field that satisfies the 
condition RabQ"^^ = 0, (A,B — z,(p), in an otherwise flat spacetime background is studied. In 
particular, it is found that the gravitational collapse of a purely gravitational wave or a null dust 
fluid cannot be realized in a flat spacetime background. The studies are further specified to the 
collapse of purely gravitational waves and the general conditions for such collapse are found. It is 
shown that after the waves arrive at the axis, in general, part of them is reflected to spacelike infinity 
along the future light cone, and part of it is focused to form spacetime singularities on the symmetry 
axis. The cases where the collapse does not result in the formation of spacetime singularities are 
also identified. 
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I. INTRODUCTION 

Gravitational collapse of a realistic body has been 
one of the most thorny and important problems in Ein- 
stein's theory of General Relativity. Particularly, in 1991 
Shapiro and Teukolsky |^ studied numerically the prob- 
lem of a dust spheroid, and found that only the spheroid 
is compact enough, a black hole can be formed. Oth- 
erwise, the collapse most likely ends with a naked sin- 
gularity in violation of Penrose's cosmic censorship con- 
jecture 0. Later, by studying the collapse of a cyhn- 
drical shell that is free of both gravitational and mat- 
ter radiation and is made of counter- rotating particles, 
Apostolatos and Thorne (AT) |^ showed analytically 
that the centrifugal forces associated with an arbitrarily 
small amount of rotation, by themselves, without the aid 
of any pressure, can halt the collapse at some non-zero, 
minimum radius, and the shell will then oscillate until it 
settles down at some final, finite radius, whereby a space- 
time singularity is prevented from forming on the sym- 
metry axis. Soon after AT's work, Shapiro and Teukolsky 
studied numerically the gravitational collapse of rotat- 
ing spheroids, and found that the rotation indeed signifi- 
cantly modifies the evolution when it is sufficiently large. 
But, for small enough angular momentum, their simula- 
tions showed that spindle singularities appeared to arise 
without apparent horizons, too. Hence, it is still possi- 
ble that even spheroids with some angular momentum 



naked singularities may be formed. However, it should 
be noted that the possibility of the formation of a horizon 
at late times cannot be excluded in the ST simulations, 
since their numerical calculations always terminate right 
at the formation of the spacetime singularities. 

Quite recently, we have studied analytically the prob- 
lem of a collapsing cylindrical shell that radiates gravi- 
tational waves and massless particles when it is collaps- 
ing, and found two physically distinguishable final states 
In one case, after the shell collapses to a minimal 
non-zero radius, it starts to expand, that is, the angular 
momentum of the dust particles is strong enough to halt 
the collapse, so that a spacetime singularity is prevented 
from forming on the symmetry axis. However, in the 
other case the rotation is not strong enough to halt the 
collapse at a finite non-zero radius, and, as a result, a 
spacetime singularity is finally formed on the symmetry 
axis. 

As shown numerically in , during the late stages of 
the collapse, the spheroid tends to fall at nearly the speed 
of light. Therefore, it is reasonable to consider the col- 
lapsing body as consisting of null fluid at its late stages 
of collapse. In this paper, we shall generalize our previ- 
ous studies 1^ to the gravitational collapse of a cylinder 
with a light-like wave- front. Specifically, the paper is or- 
ganized as follows: In Sec. II, we consider the general 
matching of two cylindrically symmetric regions across a 
null hypersurface, using the Barrabes-Israel (BI) method 
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1^, with the possibihty that on the null hypersurface a 
null dust shell may appear. It should be noted that the 
same problem was also studied in ||^ but in double null 
coordinates. As we shall see below, the method used 
there is much simpler than the BI one to be used in this 
paper, but the advantage of the latter is that it allows 
the independent choice of coordinates in both sides of 
the hypersurface, and that the use of the radial coor- 
dinate in the BI method makes physical interpretations 
of the problem more transparent. In Sec. Ill we ap- 
ply the general formulae to the case where inside the 
null hypersurface the spacetime is Minkowskian, while 
outside the shell the spacetime satisfies the condition 
RabQ^^ = 0, {A,B = In Sec. IV, we apply 

the general formulae further to the collapse of cylindri- 
cally symmetric purely gravitational waves, while in Sec. 
V our main conclusions are presented. 



symmetry be well defined, the spacetime in this region 
has to satisfy several conditions pi]-[l^: 

(i) The existence of an axially symmetric axis: The 
spacetime that has an axially symmetric axis is assured 
by the condition. 



(2.3) 



as r ^ 0+. 

(ii) The elementary flatness on the axis: This condition 
requires that the spacetime be locally fiat on the axis, 
which in the present case can be expressed as [Oj 



as r — > 0+, where X is given hy X — \\d^\\ — Igipipl^ 
Since the trajectory of the shell is null, we have 



(2.4) 



II. DYNAMICS OF CYLINDRICAL NULL 
SHELLS WITHOUT ROTATION 

Both static and dynamic |^,^ cylindrically symmet- 
ric and timelike thin shells with zero total angular mo- 
mentum have been studied previously. In this section, we 
shall give a general treatment for cylindrically symmetric 
null shells, which connects two cylindrical, otherwise, ar- 
bitrary regions. In some cases the null shell can be made 
disappear, and the corresponding null hypersurface be- 
comes a boundary surface In the following we shall 
consider the latter as a particular case of the former and 
study them all together. 

Assume that the shell, located on the hypersurface E, 
divides the whole spacetime into two regions, V^. Let 
denote the region outside the shell, and V" denote 
the region inside the shell. In V^, the metric takes the 
form 



f+{T, R)dT^ - g+{T, R)dR^ - h+{T, R)dz^ 
-l+{T,R)dcp\ {R>Ro{T)), (2.1) 



where = {T, R, z, ip} with -oo < T,z < 

+00, < ip < 27r, and R — Ro{T) is the location of 
the shell in the coordinates If the sources are con- 

fined within a finite region in the radial direction, then 
the spacetime should be asymptotically flat in the radial 
direction as 7?. ^ oo, where TZ denotes the proper radial 
distance. 

In V~, the metric takes the form 

ds'i = f^it, r)dt^ - {t, r)dr'^ - {t, r)dz^ 

-l-{t,r)d^\ (r<r„{t)), (2.2) 

where {x~^} = {t, r, z, ip} are the usual cylindrical 
coordinates chosen in the region . In this region we 
shall choose the radial coordinate r such that r = rep- 
resents the symmetry axis. The hypersurface S in this 
region is described by r = ^^{t). To have the cylindrical 



R'.{T) = e[^— 



r'oit) 



1/2 



1/2 



R=Ro(T) 



=ro{t) 



(2.5) 



where a prime denotes the ordinary differentiation with 
respect to the indicated argument, and e = ±1, with 
the sign corresponding to expanding shells, and the 
sign "— " to collapsing shells. On the shell, the intrinsic 
coordinates will be chosen as {^°} = {w, z, ip}, (a = 
1, 2, 3), where w denotes a null coordinate, which will be 
defined explicitly below. In terms of the metric on 
the shell takes the form 

ds^\^^ 5* abd^^'d^'' = ~h{w)dz'^ - l{w)dp^, (2.6) 

where ab denotes the induced metric on the hypersur- 
face. The first junction condition that the metr ics i n both 
sides of the shell reduce to the same metric (2.6) on S 
requires that 



hiw) = h+iT,Ro{T)) = h-it,roit)), 
liw) = l+{T,Ro{T)) ^ l-it,roit)), 



(2.7) 



where on the hypersurface S we have T = T{w) and 
t = t(w). The normal vector to the hypersurface E is 
given in and , respectively, by 



nt = {-R'oiT), 1, 0, 0}, 
K = {-r'oit), 1, 0, 0}. 



Following Barrabes and Israel (BI) 
the vectors Nj^ via the relations. 



(2., 



let us introduce 



0, 7V±A„± ^ 1^ (2.9) 
we find that in the present case they take the form. 
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(2.10) 



Defining tlie quantities, e^'^ by e^'^ = dx^'^/dS^'^, we 
obtain 



dT 



(1) dw 



(2) 



6(3) - Otpj 



^(2) 



dt 
dw 
5>t. 



'(3) 



5^' 



(2.11) 



To make sure that A'^,^, defined on the two faces of the 
hypersurface S, do represent the same vector TV^, we 
impose the conditions 



A (a) 



-A 

A ^(a) 



(a=l, 2, 3). (2.12) 



In the present case, it can be shown that the above con- 
ditions reduce to 



_u/2 dt 
= d^ 



(2.13) 



The above equation can be also considered as defining 
the relation between T and t on the hypersurface. Fol- 
lowing BI 1^, we define the extrinsic curvature of the 
hypersurface S as. 



A^Ae^.){<).. + r^M<)}' (2-14) 



which in the present case has the following non- vanishing 
components. 



,1/2 



d^T 1 fdrV ( /,T .9t 



dw"^ 



2\dw 




n 





fdT 


+( 


\dw 






w 








w 


!•( 


.T = 





(2.15) 



from the above expressions by the replacement 

h+, l+, Ro{T), T, R 

g^, h^, l^, ro{t), t, r. 



(2.16) 



Introducing the null vector Z° along the direction dw 
by = 5°, it can be shown that the surface energy- 
momentum tensor Tab is given by B, 



ab 



1 

2k 



r ac ibjd I bd ja jc 

[g* I I + g* i I 

J* I' I' - 5* ' ' I led, 



where k = 8ttG/c^, 

lab = 

and 



T^ab) 



(2.17) 



(2.18) 
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[g* ab] = 





-h[w) 











.0 





-l{w 


. 






■0 










[9f] = 





^h-^{w) 









.0 





-l 





(2.19) 



Inserting Eqs.(2.13), ( ^.15 ) and the corresponding ex- 
pressions for 7?.^j, into the above equations, we find that 
t"'' can be written in the form 

t"'' = pri^ + p {z'^z^ + (^''/) , (a, 6 = w, z, if), (2.20) 



where 



P = 



1 

2k 



h{w) l{w) 

1/2 



P = 




and 



Za^h'^^w)Sf„ Va = h^'^(w)5t. 



(2.21) 



(2.22) 



Cle arly, the surface energy-momentum tensor given by 
Eq.( ^.20 ) represents a Type II null fluid in the classifi- 
cations of Hawking and Ellis |jl^ , where the null fluid is 
moving along the lines defined by the tangential vector 
I"' with its isotropic pressure p in the dz and dip direc- 
tions. It is interesting to note that when /+ — 5+ and 
= g~ , we have p — Q. 
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III. GRAVITATIONAL COLLAPSE OF MATTER 
FIELDS SATISFYING RabG^^ = 

As an application of the formulae developed above, let 
us consider some specific models that represent gravita- 
tional collapse. Therefore, from now on we shall consider 
only the case where e = — 1. The wave front of a collaps- 
ing cylinder is given by R — Ro{T) in the a;"*"^ coordinates 
and r = ro(t) in the coordinates. Before the wave 
arrives, we assume that the spacetime is flat, that is. 



dt^ - dr^ - dz^ - r'^dip^, (r < ro(i)) ■ (3.1) 



Clearly, in the present case the cylindrical symmetry of 
the spacetime is well defined, and the axis r = is free 
of any kind of spacetime singularities [cf. Fig. 1(a)]. 



(a) 




(b) 



tA 








FIG. 1. The matching of a Minkowskian region M4 with an 
exterior V^: (a) In the z = Const, plane. The hypersurface 
E is a null surface represented by the circle r = ro(i). (b) 
In the {t, r) -plane. The hypersurface E now is represented 
by the straight line t + r = 0. In the region t + r < 0, the 
spacetime is Minkowskian, while in the region t + r > 0, the 
spacetime is curved. 



Then, the condition ( |2.5| ) yields, 

ro{t) - -{t-to), 



(3.2) 



where to is a constant. By a translation transformation 
in t, we can always set to = 0. In the following, we shall 
always assume that this is the case. Since ro{t) > 0, we 
must have < < on the hypersurface E. Then, we can see 
that the history of this hypersurface in the (t, r)-plane is 
represented by the part of the light cone t + r = 0, t < 0, 
as shown by Fig. 1(b). 

In the region R > Ro{T), where collapsing cylinder 
is present, the metric, without loss of generality, can be 
cost in the form 



dsl 



(3.3) 



wher e M , U an d V a re fu nctions of T and R only. Then, 
Eqs.(U), 1^ and (|2T|) yield 



RoiT) = -T > 0, V (T, RoiT)) = U (T, i?o(T)) , 



roit) 



-uiTMoiT))^ dt = e-^"^^'''<><^^^^dT. (3.4) 



Inserting the above equations into Eq.( 2.21 ), we find that 
1 



2k 



{U,T + U,R + e^}\ 



p^O. (3.5) 



To further study the problems, in the following we shall 
consider only the case where the cylinder is made of mat- 
ter that satisfies the condition 



(3.6) 



where denotes the Ricci tensor. It can be shown that 
several matter fields, such as, null dust fluid, massless 
scalar field, and electromagnetic field, sati sfy the above 
condition |l|,|l^ . For the metric given by , Eq. ( ^ ) 
yields 



U 



R 



0, 



U^TT — Urp — U^RR 4 

which has the general solution 

U = -\n [a{T + R) + b{T - R)] , 



(3.7) 



(3.8) 



where a and b are arbitrary functions of their indicated 
arguments, subject to 



a + b>Q. 

Inserting the above expressions into Eq. (| 
2a' - 1 



2K.{a + 6) 



p = 0. 



(3.9) 

]) , we find that 
(3.10) 



R=-T 



According to 



a) a'b' < 0, b) a'b' > 0, 

c) a' = 0, b' ^ 0, 

d) a' 7^ 0, 5' = 0, 



(3.11) 
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the solutions will have physically different interpreta- 
tions, thus, in the following let us consider them sepa- 
rately. 

Case a) a'h' < 0: In this case the normal vector to 



the hypersurfaces <I> 



Const, is spacelike. As a 



matter of fact, it is easy to show that 

$,a«',/35"'' = 4a'6'e*^ < 0. (3.12) 
Without loss of generality, now we can choose 

(3.13) 



a^liT + R), b = -\{T-R). 



In fact, introducing two new coordinates T and R via the 
relations, 

i? = a + 6, r = a-6, (3.14) 
we find that the corresponding metric takes the form, 
dR^) 



ds\ = e-^' {df^ - d^) - R {e'^dz^ + e-^'dip^^ 



where e' 



-M 



(3.15) 

e~^/(-4a'5'), and M and V now are func- 
tions of T and R only. Clearly, th e ab ove form of metric 
is the same as that given by Eqs.(3.3) with the choice of 
Eq.( |3.13 ), from which we find that 



e^^ = a + b = R. 



(3.16) 



Combining Eqs.(3.2) and (3.4), we find that the first 
junction condition now requires 

ro(t) = i?o(T) = -i = -r>0, 
V{T,Ro{T)) =-ln[RoiT)], Af (T, i?o(T)) = 0. (3.17) 



Substituting Eq. (|3.13|) into Eq. (|3.10|) , on the other hand, 
we find that 



p = = p. 



(3.18) 



That is, in this case the matching across the hypersurface 
E is always smooth, and no matter shell appears on the 
hypersurface. Using the C-energy defined by Thorne p7[. 



R 



(3.19) 



we find from Eq.(3.17) that on the hypersurface S we 
have 



Ce (T, i?o(r)) = o. 



(3.20) 



That is, the total C-energy inside the cylinder R = Ro{T) 
is zero. This is what one would expect, as the hypersur- 
face E now is free of matter shells and the spacetime, 
before the wave front arrives, is fiat. 

Case b) a'b' > 0: In this case it can be shown that the 



normal vector to the hypersurfaces $ = e 



Const, is 



timelike. Following a similar argument as that given in 
the last case, one can show that now we can choose 



a=-liT 



R), b^-^{T-R), 



so that 



^ a + b^ -T >0. 



Inserting Eq.(^.2l[) into Eq.(|3.10|), we find 
1 



P = 



'n{-T) 



, p = 0. 



(3.21) 



(3.22) 



(3.23) 



Since on the hypersurface E we have T < 0, from the 
above we can see that the energy density p is always 
negative. This is not physical and we shall not consider 
this case further in the following discussions. 

Case c) a' = 0, h' ^ 0: In this case, without loss of 
generality, we can set a — Q. Then, it can be shown that 
the hypersurfaces $ = — b{T — R) = Const, now 
become null, and the function b is positive and otherwise 
arbitrary. Hence, from Eq.( 3.ld| ) we find 



1 



2Kb{2Ry 



(3.24) 



from which we can see that the energy density p is always 
negative, too, and this case should be also discarded in 
the following discussions. 

Case d) a' 7^ 0, b' = 0: Similar to the last case, it can 
be shown that now the hypersurfaces $ = = Const. 
are also null, and the function a is arbitrary, subject to 
a > 0, with the choice 6 = 0. Then, the corresponding 
surface energy density is given by 



2a'(0) - 1 
2Ka(0) 



On the other hand, from Eq.( p.4[ ) we find that 



-U(T,R.o(T)) 



a(0). 



(3.25) 



(3.26) 



which is not consistent with Eq.(3.2). Thus, this case 
has to be discarded as representing gravitational collapse, 
too. 

It is interesting to note that in the last two cases the 
only possible matter field is a null dust fluid with the 
energy-momentum tensor (EMT) being given by 



Tfj^i, — pkfiky 



(3.27) 



where is a null vector. In the last case it is given by 
kf^ — Sj^ + 5^, and the corresponding EMT represents 
an ingoing null dust fluid moving along the hypersur- 
faces T + R = Const., while in Case c), it is given by 
kf^ — Sj^ — Sj}, and the corresponding EMT represents an 
outgoing null dust fluid moving along the hypersurfaces 
T — R = Const. In both of the two cases, an, ingoing 
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in the last case, and outgoing in Case c), purely gravi- 
tational cylindrical wave is present, except for the case 
where V = Const. 

In review of all the above, we find that the gravitational 
collapse in the last three cases in an otherwise flat and 
cylindrical spacetime background cannot be realized. The 
only case that is relevant to gravitational collapse is Case 
a). This result is a little bit of surprise, and mainly due to 
the requirement that, before the wave front arrives, the 
spacetime be flat and take the cylindrical form Eq.( ^.l[ ). 
Although Birkhoff 's theorem for the spherically symmet- 
ric case is not applicable to the cylindrically symmetric 
case , from the physical point of view this requirement 
seems quite reasonable. 

To further study the problem, let us now turn to con- 
sider gravitational collapse of cylindrical purely gravita- 
tional waves. 



IV. GRAVITATIONAL COLLAPSE OF PURELY 
GRAVITATIONAL WAVES 

As shown in the last section, the only case that is 
relevant to gravitational collapse in an otherwise flat 
background is Case a), in which the null hypersurface 
r = ro{t) [or R = Ro{T)] is free of matter. Tab = 0, and 
the function U{T, R) is chosen as 



U = - InR. 



(4.1) 



For this gauge, it can be shown that only three of the 
vacuum Einstein field equations Rp^^ = are independent 
and can be written in the form 



n^T = rvtVb., 
1 



V 



R ' 



V 



,RR 



R 



:V 



,R 



Vtt = 0, 



where 



n{T,R) = - \nR~ M{T,R). 



(4.2) 
(4.3) 
(4.4) 

(4.5) 



Thus, once the solution of Eq.(4.4) is known, the fun ction 
Q. ca n be obtained by the quadrature from Eqs.(4.2) and 
(^. The general solution of Eq.(^) is given by@jl9| 



V 



[A{uj) sin(wT) -I- B{uj) cos(wr)] Jq {LoR)duj 
[C{uj) sin(wT) + D{uj) cos(wT)] iVo (u;R)duj 



N 



aT + l3\nR + Y,hk In(^), 



k=l 



where 



R' 



1/2 



(4.6) 



(4.7) 



a, /3, A{u>), B{llj), C{uj), D{uj), uj, hk and Tk are ar- 
bitrary constants, and N is an integer. The functions 
Jq [ojR) and iVo {ojR) denote the Bessel and Neumann 
functions of zero order, respectively. The integration con- 
stants Ti can be real or complex, and in general takes the 
form 



Tk + i Wk, 



(fc=l, 2, 



N), 



(4.8) 



where T^'s and Wk's are real. If one starts with a complex 
term, say, /ij, then its complex conjugate /i* should be 
also included, in order to make the metric coefficients to 
be real. Note that when Tk is real, the solution is defined 
only in the regions {Tk - T)^ > i?^ [cf. Fig. 2]. 




FIG. 2. The function /ifc appearing in the expression (t.f:) 
is defined only in the regions (r^ —T)^>R^ when Tk is real. 
On the light cone, we have fik ~ R. 



On the hght cone {Tk-TY = i?^, we have ln(//fc/i?) = 
0. Then, from Eq.(4^) we can see that V{T, R) matches 
continuously to the solution that does not include /ifc in 
the region {Tk — T)^ < R? . Belinskii and Francaviglia in- 
terpreted this as gravitational shock- wave solutions , 
but later Gleiser showed that in general the light cone 
also support a null dust fluid Yet, Curir et al. found 
that sometimes the solutions have curvature singularities 
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on the light cone |22|. In any case, since here we are 
mainly concerned with gravitational collapse of gravita- 
tional waves only, we shall not consider this possibility 
any more, and simply assume that all the /i^'s are com- 
plex, as in the latter case there are no such restrictions 
and the solutions are valid in the whole spacetime. The 
term related to /ife is usually referred to as generalized 
soliton solutions [ p3| . 

On the hypersur face E where T ~ — i?, the first junc- 
tion condition Eq.(3.I7) requires 



V[-R,R) = - XnR, M{-R,R) = 0. 



(4.9) 



Since for the general solution of Eq.( ]4.6| ), the function 
M(T, R) is not known, in the following let us consider 
the above condition only in the two limits: R ^ and 
R ^ oo. When i? ^ 0, we have 



MtoR) 
NoiujR) 



In 



1, 
2 

TT 

2Tfe, 
R^ 
2Tk 



LUR 



(4.10) 



We assume, without loss of generality, that in Eq.(|4.i 
the first n ^k's take the values of and the rest {N — n) 
/ife's take the values of fi^ . Then we find that 



as R 



V{-R, R) 
0, where 



iP + h- + d) In R, 



h ^Y.^k^ d = - 



fe=i 



D{uj)duj. 



The junction condition Eq.(4.11) requires 

P + h- + d = -1. 
When i? ^ 00, we have 

2 



(4.11) 



(4.12) 



(4.13) 



JaituR) 



No{loR) 



TTLOR 
2 

TTLoR 



1/2 



COS LjR 

4 



1/2 



sin I ujR 

4 



R, 



from which wc find that 

V{-R,R) -aR + pinR. 



(4.14) 



(4.15) 



Thus, the first junction condition Eq.(4.9) in this limit 
requires that 



a = 0, (3= -I. 



(4.16) 



Combining Eqs.(4.13) and (4.16) wc find that 



0, /3: 



-1, d + h~ 



0. 



(4.17) 



To further study the problem, now let us turn to study 
the spacetime outside the hy per surface S, i.e., region 
. In particular, following [p3| , we shall concentrate 
ourselves in the four different regions: (i) The space- 
like infinity, \T\ <C R — > oo; (ii) the future null in- 
finity, R ^ T oo; (iii) the timelike future infinity, 
i? <C T ^ oo; and (iv) the focusing region, R T ^ 0^ 
[cf. Fig. 3]. 

(i) The spacelike infinity, \T\ <^ R oo: In this limit, 
let us first consider the first two terms in E q.([1.6| ) , which 
will be denoted by Vp b {T, R) . Using Eq. ( ^1.14| ), we find 
that I2I 



Vfb 



1 



2ttR 



1/2 



{E+{uj) cos [uj{R -T) + (j)+{Lj)] 



duj 



+E-iLo) cos [^(i? + T) + rito)] } (4.18) 
where 

,^ _ fr ,^ _i_ n,', ,m2 , r/^/, .\ ^ d/, ,m2-i V2 



E^iio) = {[A{u;) ± Dico)]' + [C(c.) T 5(^)1'} 



(f>^{uj) = tan 



-1 



C{u;)tB{uj) 



A{uj)±D{uj) 



(4.19) 



It is clear that, as R ~f 00, these two terms represent a su- 
perposition of two traveling gravitational waves along the 
radial direction, one is ingoing represented by the term 
E~ {uj) and the other is outgoing represented by the term 
E'^{uj). When the waves have only a single wavelength, 
say, Ao, it can be shown that the corresponding function 
M{T,R) takes the form [||, 

Mfb -^^lnR-^[{Al + Bl + + DI) R 
+ 2{BoCo - AoDo)T] 
- i- {Ea cos [2a>o(r + R)- 6^] 



+ E+cos[2ujq{R^T)-9+]} 
where ujq = 2tt/\q, 



9± 
^0 



tan 



2(Ao±i5o)(CoT5o) 



(Ao ± DoY - (Co T B^f 



(4.20) 



(4.21) 



and Aq = A{u)q), etc. It can be shown that the space- 
time is singular in this limit. Since we are considering 
gravitational collapse, the more interesting case would 
be that the spacetime is asymptotically flat in the radial 
direction, and it is possibly singular only at i? = for 
T > 0. Then, we are sure that such singularities are in- 
deed formed by the collapse. Therefore, in the following 
we shall assume that these two terms vanish identically, 



A{uj) = B{uj) = C{lu) = D{lo) = 0. 



(4.22) 



Hence, fro m Eq .( [I.12| ) we find that d = 0. Substituting 
it into Eq.(4.17) we obtain 
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h- = j2h^ = o. 



(4.23) 



fc=i 



The third term in Eq.(4.6) vanishes too, because of 
the first junction condition (4.17). Thus, in the fohowing 



let us concentrate ourselves only on the fourth and last 
(soliton) terms with /3 — —1. When the soliton term 
vanishes, the spacetime becomes Minkowskian. 

As we mentioned previously, in this paper we consider 
only the case where Tk is complex. For the complex Tk it 
is found convenient to introduce the following quantities 
p3|, 



where 



Uk 



Lk- 

i?2 



A^fe = — = R<^k 



(4.24) 



1)1/2, 




22.0 




i?2 + 




Tk = Tfe" 





1/2 



lk 



R{1 + (Tk) 

If we further assume that 
h 

h 



(4.25) 



fc+n/2 ^\^k) =^k -iPk^ 



k+{N+n)/2 ~ \"'k+nj ~ "fc ^ h^k i 

where a^'s and /3^'s are real constants, we find that 

N 

\nR + J2hkln[^) 



(4.26) 



soliton 



Ini? 



where s = N/2, and 



fc=i 



n 

fc=i 



gk 



In (Tfc + /fc7fc 



(4.27) 



9k 



-^o^k^ 

/c+n/2 ' 



fk 



1 < fc < n/2 

n/2 + 1 < fc < 7V/2 ' 

1 < fc < n/2 

n/2 + 1 < fc < iV/2 ■ 



(4.28) 



In terms of a^, 's, Eq.( 4.23| ) takes the form 



n/2 

fc=i 



(4.29) 



The general solution (4.27) was studied in details in 
psf , thus in the following wc shall only summarize its 
main properties, and for the details, we refer readers to 
p3[ . In particular, it can be shown that 



<Jk^l+^-^+0{R-% 
7.-f-| + 0(i?-^), 



(4.30) 



as i? — > cxD. Thus, from the above we can see that the 
terms that are proportional to ak are negligible, while 
the terms that proportional to 7fc have the following dis- 
tributions to the metric coefficients. 



where 



fe=i 



^soliton 
soliton 



fk 



Ini? + ^ fklk, 



fc=l 



-UnR, 




{N+n)/2 

E PI 

fe=(n+2)/2 



(4.31) 



(4.32) 



Thus, in this limit the solution is singular unless / = 
With the same reason as that given following Eq.( 4.21 ) 
in the following we shall assume that 



ri/2 

fc=i 



(N+n)/2 

E Pt 

fe=(n+2)/2 



0, (/ = 0). 



(4.33) 



Then, we can see that the solution will represent the 
gravitational collapse of purely gravitational cylindrical 
waves with lightlike wave-front, on which no dust shell 
appears. The spacetime is asymptotically flat in the ra- 
dial direction and locally flat near the symmetry axis. 

(ii) The future null infinity, R ^ T ^ oo: In this limit, 
it can be shown that 



CTfc 

lk 



1 + O \R-^I'^ 
n + 0(R-'/^ 



(4.34) 



and that the metric coefficients behave like 
V{T, R) ^ - InR + O {r-^I^ 
M{T,R) O (i?"i/2 



(4.35) 



Thus, in this region the spacetime is asymptotically 
Minkowskian, too. A detailed study of the Riemann ten- 
sor showed that the spacetime in this region is Petrov 
type N and represents an outgoing gravitational wave 
along the light cone ||2^ . 

(iii) The timelike future infinity, i? T — + oo: The 
study of the solutions in this region will give us some 
information about the final state of the collapse. It can 
be shown that 



lk 



4T2 



(4.36) 
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and 



^soliton 
soliton 



-(1 + 5+) In i? + .g+ InT, 
-ig+(2 + .g+)lni? 



where 



{N+n)/2 



-2 E 



(4.37) 



(4.38) 



fe=l+n/2 



( Hi) 



(iv) 




(i) 



FIG. 3. The corresponding Penrose diagram for gravita- 
tional collapse of purely gravitational cylindrical waves. The 
wave-front is denoted by the straight line T + R = 0. be- 
fore the waves arrives, the spacetime is flat, denoted by 
M4. The spacetime is asymptotically flat at spacelike infinite 
|T| <C -R ^ 00 in the radial direction. When the wave- front 
arrives at the symmetry axis {T, R) = (0, 0), in general part of 
the wave is reflected to inflnity along the light cone T~ R = 
and part of it is focused to spacetime singularity on the axis, 
denoted by the curved vertical line. 



Then, we find that the spacetime in this region is al- 
ways singular, except for the cases 



When g+ = 0, we can see that the collapsing gravita- 
tional waves are reflected completely from the symmetry 
axis to spacelike infinity along the future light cone, and 
nothing is left behind, so that the final state of the space- 
time is Minkowskian. 

(iv) The focusing region, R T —^ 0^: This is the 
region where the collapsing gravitational waves just ar- 
rive at the axis, and start to accumulate there. It can be 
shown that 



Ik 



(THY 



i?2 



(4.40) 



and 



' soiiton 



'-soLiton 



{l + g+)lnR, 



:.9+(2 



In R, 



(4.41) 



(i) < g+ < 1, or (ii) g^ 



(4.39) 



where 7^ is a constant. Thus, the spacetime in this region 
is singular, unless the conditions given by Eq.(4.39) hold. 

In review of all the above, we have the following: 
The only solutions that represent gravitational collapse 
of purely gravitational cylindrical waves in an otherwise 
flat spacetime background are g iven by Eg. ( 4. 27 ) sub- 
jected to the conditions of Eqs.(O^) and ( 4.33 ). Af- 
ter the waves arrive at the symmetry axis, part of them 
is reflected to spacelike infinity along the light cone 
i? ~ r — > 00, and part of them is focused to form 
spacetime sin gula rity on the axis, except for the cases de- 
fined by Eq.(4.3!;). The corresponding Penrose diagram 
is given by Fig. 3. 



V. CONCLUSIONS 

In this paper, the general matching of two cylindri- 
cally symmetric regions across a null hypersurface is 
studied using the Barrabes-Israel method. It is shown 
that in general the hypersurface supports a null dust 
shell. As an application of these formulae, we study 
the collapse of a matter field that satisfies the condition 
RabQ^^ = 0, {A,B = z,ip), in an otherwise flat space- 
time background. It is found that the only case that 
represents gravitational collapse is that where the nor- 
mal vector to the hypersurface <I> = Const, is space-like, 
where $ = {gzzgipipY^^ ■ As a result, the gravitational 
collapse of a purely gravitational cylindrical wave or a 
pure null dust fluid cannot be realized in a flat spacetime 
background. 

The general formulae are further applied to the col- 
lapse of purely gravitational cylindrical waves in a flat 
spacetime background and the general conditions for such 
collapse are found. It is shown that after the waves arrive 
at the axis, part of them is reflected to spacelikc infinity 
along the future light cone, and part of it is focused to 
form spacetime singularities on the symmetry axis. The 
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only cases where n o spa cetime singularity is formed are [23] J. Garriga and E. Verdaguer, Phys. Rev. D36, 2250 
these given by Eq. (|4.39D . (1987) . 

[24] P.J. Adams et al, Astrophys. J. 253, 1 (1982). 
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